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The conductivity of graphene on a boron nitride substrate exhibits features in the terahertz (THz) 
and infrared (IR) frequency regimes that are associated with the periodic moire pattern formed by 
the weakly coupled two-dimensional materials. The THz and IR features are strongest when the 
two honeycomb lattices are orientationally aligned, and in this case are Pauli blocked unless the 
Fermi level is close to ±150 meV relative to the graphene sheet Dirac point. Because the transition 
energies between moire bands formed above the Dirac point are small, ac conductivity features in 
n-doped graphene tend to be overwhelmed by the Drude peak. The substrate-induced band splitting 
is larger at energies below the Dirac point, however, and can however lead to sharp features at THz 
and IR frequencies in p-doped graphene. In this Letter we focus on the strongest few THz and IR 
features, explaining how they arise from critical points in the moire-band joint density-of-states, 
and commenting on the interval of Fermi energy over which they are active. 


Introduction :— Strong interactions between graphene 
sheets and light have motivated research aimed at poten¬ 
tial graphene-based optoelectronic deviceJ^lf^ operating 
in the terahertz (THz) and the infrared (IR) frequency 
regimes, a range of the electromagnetic spectrum that 
is important for imaging, sensing, and communications 
technologies.lS®^ The ability to tune carrier density in 
graphene using gate voltagei^^ is a great advantage in 
most application possibilities. 

The ac conductivity of a neutral graphene sheetpHUl is 
nearly frequency-independent (Tq = 7re^/2/i; adding n or 
p type carriers transfers oscillator strength over the fre¬ 
quency interval (0, 2u)f) which is Pauli blocked for inter¬ 
band transitions to a Drude peak of equal weight. When 
a graphene sheet is aligned with a hexagonal boron ni¬ 
tride substrate secondary Dirac points are induced by the 
substrate and, because they are gapped, lead to a large 
joint density-of-states for low-frequency optically active 
transitions. Like the interband conductivity of an iso¬ 
lated graphene sheet, the sharp conductivity features as¬ 
sociated with these transitions can be turned off and on 
by adjusting the position of the Fermi level. (See Fig. [^. 
In this Letter we present a theory of the conductivity at 
Fermi level values close to secondary Dirac points. 

Hexagonal boron nitride is a wide band gap 
semiconductoi^i^ with a hexagonal lattice structure and 
weakly coupled layers. Placing graphene on a hexago¬ 
nal boron nitride substrate retains the high quality of 
graphen^isE^ while modifying its band structure and 
therefore its optical response. The substrate’s pattern of 
negatively charged boron atoms and positively charged 
nitrogen atoms on opposite honeycomb sublattices alters 
the graphene 7r-band Hamiltonian. If the two honeycomb 
lattices had perfectly matched lattice constants and ori¬ 
entations, the primary effect of this intera ction would be 
to simply to open a gap at the Dirac point .1121^ However, 
the boron nitride lattice constant is about two percent 
larger than that of graphene and differences in orien- 
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Figure 1. Pauli blocking in bare graphene, (a) This car¬ 
toon shows the Dirac cone of doped isolated graphene with a 
Fermi level ef given by the dashed line. Blue arrows represent 
interband transitions from the valence band to the conduc¬ 
tion band. The middle arrow, with frequency 2ujf = 2eF/h, 
is at the minimum frequency for interband transitions. Be¬ 
low 2ujf, transitions are not allowed due to the Pauli exclu¬ 
sion principle, (b) Density of states of aligned graphene on 
hBN. The shaded regions represent the locations of the sub¬ 
strate induced gapped Dirac cone replicas. The conductivity 
is strongly Fermi level dependent near these energies because 
of Pauli blocking. 

tation are typical in exfoliated samples,leading to 
a more complex electronic struct ure. In nearly aligned 
layers, long-period moire pattern^ ^^bspz i form and influ¬ 
ence all physical properties. Experiments show that a 
gap opens at the Fermi level of neutral graphene sheets, 
with a value that is dependent on electron-electron in¬ 
teractions, and on strains induced by the lattice con¬ 
stant and orientation mismatch .I22H2I1 At perfect align¬ 
ment (zero relative rotation angle between the graphene 
and boron nitride), the moire wavelength is around 15 
nm. In addition to the gap at ch arge neu trality, sec¬ 
ondary gapped Dirac points appeai|2IE2E6H2Sl ^t Fermi 
energies corresponding to ±4 electrons per moire period, 
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Figure 2. Band structure and ac conductivity of graphene on boron nitride, (a) Band structure of orientationally 
aligned graphene on boron nitride. The dashed lines show the Fermi levels for which the ac conductivity has been calculated, 
(b) Schematic picture of the moire Brillouin zone, outlined in red. The bands are plotted along the black lines in momentum 
space. Blue dots represent the first shell of moiree reciprocal lattice vectors, (c) The ac conductivity when the Fermi level is 
at the Dirac point ef = 0. The red dashed line shows that (j(tu) approaches the bare graphene value cro = 7re^/2/i at high 
frequencies, (d) The ac conductivity when the Fermi level is in the conduction band and near the secondary Dirac point with 
Ef = hvG/'Jh (black) and near the secondary Dirac point with ef = hvGl2 (blue). G is the magnitude of the moire reciprocal 
lattice vectors in the first shell, (e) The ac conductivity when the Fermi level is in the valence band and near the secondary 
Dirac point with ef = —hvG/VS (black) and ef = —hvG/2 (blue). 


which corresponds in perfectly aligned graphene to Fermi 
levels ^ ±150 meV. This Fermi level scale is set by the 
moire period, and can easily be reached by gating. The 
substrate interaction yields particle-hole asymmetry, and 
in particular leads to very diff erent densiti es-of-states in 
conduction and valence bands 

The ac conductivity is dominated by the intraband 
Drude peak and interband features associated with sin¬ 
gularities in the joint density-of-states. It is therefore 
strongly affected by the band structure modifications in¬ 
duced by the boron nitride substrate. Previous work has 
highlighted the feasibility of using optical absor ption to 
determine properties of the substrate interaction.l22l311 In 
this paper, we use a substrate interaction HamiltoniarP^ 
derived from ab initio electronic structure calculations 
to evaluate the ac conductivity of graphene on boron ni¬ 
tride. We focus on the case of perfect orientational align¬ 
ment. In contrast to other substrates for which the ac 
conductivity in the THz regime is dominated by a broad 
intraband Drude peak, the moire pattern formed by 
graphene on boron nitride sometimes induces sharp THz 
peaks due to transitions between Bloch bands formed by 
the moire superlattice. The particle-hole asymmetry of 
the moire Bloch bands is strongly reflected in the THz 
and IR conductivity which is always Drude-dominated 
when the Fermi level lies above the Dirac point, but is 


interband-dominated when the Fermi energy lies in a rel¬ 
atively narrow interval below the Dirac point. The qual¬ 
itative change in THz and IR conductivity with Fermi 
energy suggests a potential mechanism for electrically 
tunable optical properties. 

THz conductivity calculation :— We employ the moire 
band Hamiltonian for graphene 7r-band electrons de¬ 
scribed in Ref IdOl in which a local periodic substrate in¬ 
teraction term is added to the k • p continuum Dirac 
model of an isolated graphene sheet. The substrate in¬ 
teraction is extra cted from ab initio electronic structure 
calculation^^nOni and is accurate for moire patterns with 
spatial periods much larger than the graphene sheet lat¬ 
tice constant.^ In previous work, we have used this ap¬ 
proach to show how strains in the graphene lattice and 
the substrate together with electron-electron interactions 
control the size of the gap at the Fermi level of neutral 
graphene sheets.!^ We have also explained the subtle way 
in which substrate induced changes in carbon-site ener¬ 
gies and hopping strengths combine to yield surprisingly 
strong particle-hole asymmetry that is manifested both 
in the density-of-states and in dc transport properties 

The real part of the ac conductivity can be evaluated 
at zero tempera ture and f requency lu using the Kubo for¬ 
mula expressio]J2I^2inni37| 
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where r is taken to be the momentum relaxation time, 
ep is the Fermi energy, ^ is the energy of moire band 
n at wave vector k, and \n, k) is the eigenstate. The step 
functions, 0(x), ensure that transitions occur only when 
one state is filled and one is empty, while the Dirac delta 
function 6(x) enforces energy conservation. The matrix 
elements can be simplified by noting that the Hamilto¬ 
nian depends on k only through the continuum Dirac 
part, hvk ■ T, where t“ are Pauli matrices (a = x,y). 

The results of a calculation of the ac conductivity in 
which Eq. was evaluated numerically are illustrated in 
Fig.H which also plots the moire bands and summarizes 
the geometry of the moire Brillouin zone. The plotted 
bands describe states near the microscopic K Dirac point 
and given eigenenergy as a function of momentum mea¬ 
sured from that point in momentum space. The Hamil¬ 
tonian near the microscopic K' Dirac point is related to 
the Hamiltonian at the K point by time reversal sym¬ 
metry; the bands and energy eigenstates are therefore 
degenerate in the absence of a magnetic field. An ex¬ 
ternal magnetic field breaks the time reversal symmetry, 
and the ac conductivity contributions from the two val¬ 
leys will no longer be identical. In this paper, we consider 
only the case of zero magnetic field. 

Because the moire pattern is triangular, we choose to 
label symmetry points in the moire band structure plots 
conventionally with subscript to, although all indicate 
points in momentum space close to the Brillouin-zone 
corner K of the microscopic reciprocal lattice. The low- 
energy electronic structure consists of moire bands asso¬ 
ciated with each of graphene’s two valleys. The Kramers 
partner of a moire band state associated with one valley 
lies in the opposite valley. For this reason the moire band 
structure of a particular valley is not time reversal invari¬ 
ant, allowing the energies at moire Km and K^ points 
to be different. When the Fermi level is at the Dirac 
point a gap opens due to sublattice symmetry breaking. 
The conductivity is therefore zero for frequencies below 
the gap, peaks above the absorption threshold, and ap¬ 
proaches th e univ ersal conductivity of isolated graphene 
(To = 7re^/2/iSl^ at high frequencies. There is no trans¬ 
port gap at any Fermi level value in the conduction band, 
so there is always a significant contribution to the con¬ 
ductivity from the Drude term peaked at w = 0. Small 
features superimposed on the Drude peak appear at finite 
oj and are associated with the onset of interband transi¬ 
tion channels as the frequency increases. In contrast to 
the case in the conduction band, for clean enough sam¬ 
ples, there is a transport gap in the valence band when 
the hole density is four per moire peirod. When the Fermi 
level lies in or just above or just below the moire band 
edges (as shown in Figure]^ the Fermi surfaces are small, 
and so the Drude contribution remains smaller than in 


the conduction band case. Interband transitions are then 
dominant. Below we associate different peaks with tran¬ 
sitions between moire bands at different points in mo¬ 
mentum space. 

TKz conductivity analysis :— In Fig. [^we have separated 
the interband conductivity into contributions from pairs 
of moire bands which cover energy intervals below the 
Dirac point energy. For each pair we have identified 
the dominant features, and associated them with tran¬ 
sitions at or near high-symmetry points in the moire 
Brillouin zone which are marked by arrows in the band 
structure panels of Fig. First consider Fermi level 
ep = El = h.vGj^/'i at which all three interband tran¬ 
sitions are active. We order the valence bands in reverse 
order of energy so that the first valence band is the one 
closest to the Dirac point in energy. 

We see in Fig. that the second subband is flatter 
over a wider region of momentum space than either the 
first or the third valence band. Accordingly third to 
first transitions make a smaller contribution to the in¬ 
terband conductivity than either third to second or sec¬ 
ond to first transitions. The two strongest features are a 
third to second transition from near the Km-, K'm points 
with energy 0.017 eV (red in the right panel) and second 
to first transitions along the Km, Mm line with energy 
0.032 eV(purple in the left panel). All features can be 
identified with particular high symmetry points or lines 
in the moire Brillouin zone, as shown by the matching 
colored arrows in the band structure plots. Generally 
speaking the strongest features are associated with in¬ 
terband transitions near the Km and K'm points in the 
moire Brillouin zone. The top row of Figure shows 
differences in band energy AA as a function of position 
in the moire Brillouin zone as a contour plot. Up to a 
matrix element factor, the optical conductivity is pro¬ 
portional to the joint density of states. For transitions 
from the 2nd to 1st band AA has a saddle point at Km 
which leads to a divergence in the joint density of states. 
In contrast, extrema in the transition energy yield jump 
discontinuities in the joint density of statepSl and weaker 
optical conductivity features. We note that in an experi¬ 
ment, the presence of impurities (which we have ignored 
for interband transitions) will broaden the features, and 
lead to a finite conductivity at saddle points. 

The relatively small change in Fermi level to ep = 82 = 
hvG/2 leads to a very different conductivity profile. Be¬ 
cause both the 2nd and 3rd bands are full, transitions 
between these two are forbidden, leading to a zero con¬ 
tribution from this pair of bands. The dominant fea¬ 
ture becomes the 2nd to 1st band transitions at the K'm 
point. This feature increases in strength as the Fermi 
level moves into the transport gap, and will decrease to 
zero as the Fermi level is moved into the first valence 
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Figure 3. Valence band conductivity features association with transitions at high symmetry points in the moire 

Brillouin zone. The interband transition energy (first row), band structure (second row) and conductivity (third and forth 
rows) for interband transitions between three pairs of valence bands. The left column corresponds to transitions from the 
second to first valence band; the middle column corresponds to transitions from the third to first valence band; and the right 
column corresponds to transitions from the third to second valence band. In the band structures (second row) dominant 
transitions have been identified with colored arrows, and their corresponding energies are marked as vertical dashed lines on 
the conductivity plots (third and fourth rows). The two Fermi energies at which calculations were performed, marked by black 
horizontal dashed lines on the band structure plots, are £1 = —hvG/y/i and £2 = —hvGl2. The third row corresponds to Fermi 
energy ef = £ 1 , while the fourth row corresponds to Fermi energy ef = £ 2 - 


band. Note that this quite large change in the conduc¬ 
tivity is seen over a relatively small change in Fermi level 
of ~ 20 meV. 

For both Fermi levels, transitions from the 3rd to 1st 
bands are relatively weak. The reason for this is two fold. 
Firstly, the magnitude of the energy difference is larger. 
Since the conductivity is inversely proportional to this 
energy difference, this leads to a smaller conductivity. 
Secondly, as shown in the middle contour plot in Fig¬ 
ure]^ the transition energies at both Km and K'm for this 
pair of bands are minima, not saddle points, and there¬ 
fore do not lead to divergences in the conductivity. The 
contribution to conductivity due to these transitions are 
jump discontinuities, and have a relatively small magni¬ 
tude compared to the saddle point divergences discussed 
above. 

Conclusions — We find that the conductivity of graphene 
on a boron nitride substrate is extremely sensitive to 
Fermi level. The system’s particle-hole asymmetry is 


strongly manifested in the ac conductivity. For Fermi 
levels in the conduction band, the Drude peak domi¬ 
nates at low frequencies. For Fermi levels in the valence 
band, on the other hand, interband transitions dominate 
when the carrier density is close to four holes per moire 
unit cell. All features in the ac conductivity are asso¬ 
ciated with transitions at high-symmetry points in the 
moire Brillouin zone which support critical points in the 
joint density of states. Depending on the Fermi level, 
transitions between the first and second moire valence 
bands, or between the second and third bands, which 
have transition energy saddle points and associated di¬ 
vergent joint densities-of-states, dominate the total con¬ 
ductivity. These large and easily tunable changes in con¬ 
ductivity may be valuable for THz or IR applications. 
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